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Abstract 

We consider within-host virus models with n > 2 strains and allow mutation between the strains. If 
there is no mutation, a Lyapunov function establishes global stability of the steady state corresponding 
' to the fittest strain. For small perturbations this steady state persists, perhaps with small concentrations 

- of some or all other strains, depending on the connectivity of the graph describing all possible mutations. 

Moreover, using a perturbation result due to Smith and Waltman [M], we show that this steady state 
PQ \ also preserves global stability. 

u' 

d ! 1 Introduction 

' ^ ' The study of the dynamics of within-host virus disease models has been a very fruitful area of research over 

O^. the past few decades. Of particular importance has been the work on mathematical models of HIV infection 

by Perelson and coauthors [101 E] and Nowak and coauthors [3] . It has spurred more recent research by 
among others Hal Smith with one of us [T] , [TB] and [T^] . 

For single-strain virus models, the understanding of the global behavior has been largely based on the 
fact that they are competitive [T] and the use of particular mathematical tools developed by Muldowney; see 
Li and Muldowney [8] for an application of these techniques to the classical SEIR model in epidemiology. 
. Of course it is well known that for globally stable systems there is a Lyapunov function, but finding such a 

function is often difficult, as illustrated by the following quote from Smith and Waltman's classical work on 
chemostats [13] on p. 37: 



Considerable ingenuity, intuition, and perhaps luck are required to find a Liapunov function. 

One of the purposes of this paper is to find such Lyapunov functions for various within host virus models 
following the ingenuity from and Another purpose of the paper is to investigate what happens if we 
. , include mutation effects in the model by allowing different virus strains to mutate into each other. This is 

' very relevant in the context of HIV where mutations have profound impact on treatment, see for instance 

[12j where a two-strain model is considered. 

Mathematically we will treat the model with mutations as a perturbation of the original model. It turns 
out that the structural properties of the mutation matrix that describes the possible mutations (in particular, 
whether this matrix is irreducible or not), dictate which single strain steady states of the unperturbed model 
persist in the perturbed model, and which don't. An obvious problem is to determine if the globally stable 
single strain steady state of the unperturbed model persists. We will show that this is always the case, 
regardless of the mutation matrix. Moreover, taking advantage of the perturbation result developed by 
Smith and Waltman in [M], we will show that this steady state remains globally stable for small values 
of the mutation parameter. In order to apply this perturbation result we will first need to establish a 
particular persistence property, uniform in the perturbation parameter, and to achieve this we invoke the 
theory developed by Hutson [51 13] , see also , which uses the notion of an average Lyapunov function. It 
will be shown that a rather simple -in fact, linear- average Lyapunov function exists. 
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The paper is organized as follows. In Section[2]we present a Lyapunov function to establish global stability 
of the disease equilibrium of a single-strain virus model. This is extended in Section [3] to a global stability 
result for a multi-strain model which does not include mutations. In biological terms, we demonstrate that in 
the absence of mutations the fittest strain of the virus drives all other viral strains to extinction. In Section 
|4]we investigate what happens if mutations are taken into account for two different models. Finally, in the 
Appendix we extend all our results to a slightly modified model which includes an often neglected loss term 
in the virus equation. 



2 Single-strain 

In this paper, we consider the basic model of the form 

t = f{T) - kVT 
T* = kVT - (3T* 

V = Nf3T* - 7^, (1) 

where T, T*, V denote the concentrations of uninfected (healthy) and infected host cells, and free virions, 
respectively. Equations ([T]) describe a general viral infection where the viral replication is limited by the 
availability of target cells T. In this model, we assume that all infected cells T* are virus-producing cells, 
that is, we do not include any intermediate stage(s) corresponding to latently infected cells. In addition, 
we do not explicitly consider the impact of the immune response. Implicitly, the immune response can be 
accounted for by the removal term — /3T*. The rate of viral production is assumed proportional to the 
removal of infected cells. In case of lytic viruses, N represents the average burst size of a single infected cell; 
whereas in case of budding viruses, N can be thought of as the average number of virions produced over a 
lifetime of an infected cell. For different infections, the actual class of the target cells in ([T]) may vary from 
the CD4:+ T lymphocytes (in case of HIV), to the epithelial cells (in case of Influenza), to the red blood 
cells (in case of Malaria). The T, T*, V notation is adopted from the classical HIV model [TP. 

All parameters are assumed to be positive. The parameters (3 and 7 are the removal rates of the infected 
cells and virus particles respectively. Following [HI |9] , we neglect the term in the F-equation that represents 
the loss of a virus particle upon infection. But all subsequent results hold when this loss term is included, 
in which case the ^-equation reads: 

V = N(3T* -jV - kVT. 

These results will be presented in the Appendix. 

The growth rate of the uninfected cell population is given by the smooth function f{T) : K-|_ M, which 
is assumed to satisfy the following: 

3ro>0 : /(T)(r-To) <0, T / Tq. (2) 

Since continuity of / implies that /(To) = 0, it is easy to see that 

^0 = (To, 0,0), 

is an equilibrium of ([1]). Effectively, Tq is the carrying capacity for the healthy cell population. 
A second, positive equilibrium may exist if the following quantities are positive: 

rp T rp^ / (T) / (T) 

Note that this is the case if and only if fij^) > 0, or equivalently by ^ that T = < Tq. In terms of 
the basic reproduction number 

n kN Tq 

:= To = ^, 

7 To' 

existence of a positive equilibrium is therefore equivalent to 7^° > 1. We assume henceforth that 7?," > 1 and 
denote the disease steady state hy E = (T, T*,V). Let us introduce the following sector condition: 

(C) if{T)-f{T))fl-^)<0. 



2 



Note that this condition is satisfied when /(T) is a decreasing function, independently of the value of T. 
For instance, [5] considers f{T) = c\ — C2T, where q are positive constants. Another example [TT] is 
/(T) = s + rT(l - ^) provided that /(O) = s > /(T). 

Theorem 1. Let (C) /lo/rf. T/ien the equilibrium E is globally asymptotically stable for ^ with respect to 
initial conditions satisfying T*(0) + V{0) > 0. 

Proof. Consider the following function on int(M'^): 

W = I [l--]dr+ r (l~—]dT+Jl-r(l--]dT. 



T I .]rr, \ T I N (3 ./l7 \ T 



Then 



W = {f{T)^kVT)[l~^+{kVT-PT*)[l-^^+^{N(iT* -^V)[l-^ 
= /(T) (1 - y j + kVf - kVT— + (if* -f]T*--^V+ 

Since from Q we have that /3T* = fct^T = it follows that 

/ f\ - V - f*VT - - VT* V 

- f f\ / f\ - V - f*VT - - VT* V 

= {f{T) - f{T)) il^-\+pT*{l~-\+PT*^- PT*^;^ + PT* -PT*—- /3T*^ + PT* 



{f{T) - f{T)) ( 1 - ^ ) - /3T* 



T T*VT VT* 
■7^ + + 



T T*VT VT* 



The first term is non-positive by (C). The second term is non-positive as well since the geometric mean of 3 
non-negative numbers is not larger than the arithmetic mean of those numbers. Hence, Vl^ < in int(R'^), 
and the local stability of E follows. Notice that W equals zero iff both the first term and the second term 
are zero, and using (C), this happens at points where: 

T f*V 

— = 1 and ^ = 1. 

T T*V 

Then LaSalle's Invariance Principle [7] implies that all bounded solutions in int(K'^) converge to the largest 
invariant set in 

M = {(T, T*, V) G int(R3 ) | ^ = 1, |^ = !}• 

Firstly, boundedness of all solutions follows from Lemma [1] which is proved later in a more general setting. 
Secondly, it is clear that the largest invariant set in M is the singleton {E}. Finally, note that forward 
solutions starting on the boundary of Ri^ with either Ti(0) or Vi(0) positive, enter int(Mi^) instantaneously. 
This concludes the proof. □ 

3 Competitive exclusion in a multi-strain model 

Let us now consider a multi-strain model: 



f{T)-J2k^V^T (4) 



i=l 

f* = hVT-P,T*, i = l,...,n (5) 
V = N,(3,T*-j,V, i = l,...,n (6) 

where all parameters are positive. Similar calculations as in the single-strain model show there is a unique 
disease-free equilibrium Eq = (Tq, 0, 0). For each i, there is a corresponding single-strain equilibrium Ei with 
positive T, T* and Vi components and zero components otherwise if and only if 
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Here, TZ'^ is the basic reproduction number for strain i which is defined by 
The positive components of are then given by 

We assume that ah Ei exist and assume without foss of generahty (by possibly reordering components) that 

< f 2 < . . . < T"-i < T" < To, (8) 

or equivalcntly, that 

i<K< <_i <---<K<nl (9) 

and will prove the following competitive exclusion principle. It asserts that the strain with the lowest target 
cell concentration at steady state (or equivalcntly, with highest basic reproduction number) will ultimately 
dominate, provided that such strain is present initially. 

Theorem 2. Assume that all Ei exist for Q — (0), that (C) holds with instead of T , and that © 
holds. Then Ei is globally asymptotically stable for @ — ® with respect to initial conditions satisfying 
T^{0) + Vi{0) > 0. 

Proof Consider the following function on U := {{T,T^ , . . . ,T*,Vi, . . . , Vn) e M2"+^ | T, T^, Vi > 0}: 



Then 
W 



n , 



Notice that the first three terms can be simplified in a way similar as in the proof of Theorem [1] and using 
the expression for in ([7]), we find that 



W = (/(T)-/(Ti))(l-^) -/3iT; 



T 



T 



T*ViT ViT^ 



T*ViT^ ViT* 



-3 



i=2 



Each of the first two terms is non-positive as was shown in the proof of Theorem [T] The third part is also 
non-positive by ([8|). Thus < 0, establishing already stability of Ei. An application of LaSalle's Invariance 
Principle shows that all bounded solutions in U (and as before, boundedness follows from Lemma [1] which 
is proved later) converge to the largest invariant set in 



(r,r*,...,T„*,yi,...,y„)ec/| ^ = 1, 



T*Vi 
T*Vi 



1, = 0, i>2 



which is easily shown to be the singleton {Ei}. Finally, solutions on the boundary of U with ri*(0)-|-Vi(0) > 
enter U instantaneously, which concludes the proof. □ 
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4 Perturbations by mutations 



In this section we expand model fl]) — ((6]) to account for mutations between the n strains. In fact, we will 

study two different ways in which mutations occur. Our first extended model can be written compactly as 
follows: 

T = f{T)-k'VT, reM+ (10) 

f* = P{n)KVT-BT*, T* eRl (11) 

V = NBT* -TV, V (12) 

while the second is written as 

T = /(T) - k'VT, T G R+ (13) 

f* = KVT-BT*, T* eRl (14) 

V = P{^i)NBT* - TV, V eRl, (15) 

In both models K = diag(fc), B = diag(/3), N = diag(A^) and T = diag(7), and the matrix P(/i) with 
G [0,1] is defined as follows: 

Pi^I) = 1 + fiQ, 

where Q is a matrix with qij > if strain j can mutate to i (for i ^ j) so that different magnitudes of qij 
reflect the possible differences in the specific mutation rates. The diagonal entries of Q are such that each 
column of Q sums to zero. Notice that P is a stochastic matrix (all its entries are in [0, 1] and all its columns 
sum to one) provided that /i < —l/qu for all i (which is assumed henceforth), and that P(0) = /. 

Lemma 1. Both system (TQ ~ (HH and (HH — (j_?5p are dissipative, i.e. there exists a forward invariant 
compact set K C R^"^"'^ such that every solution eventually enters K . 

Proof. From ^ and ^ follows that t < f{T), hence 

limsupT(i) < To, (16) 

t — >oo 

provided solutions to both systems are defined for all t > 0. To see that this is indeed the case, we 
argue by contradiction and let {T{t), T*{t), V{t)) be a solution with bounded maximal interval of existence 
:= [0,tmax)- Then necessarily T{t) < max(T(0),ro) ■= Tmax for all t G This implies that on X+, the 
following differential inequality holds for the solution of system (fTOl) — p2|) : 



T* < Pi^i)KVT^,^-BT* (17) 

V < NBT* ~ TV, (18) 

or for system (fT3| — (fTS)) 

t* < KVT^^^-BT* (19) 

V < P{pC)NBT* - TV, (20) 

respectively. Notice that the right hand sides in the above inequalities are cooperative and linear vector fields. 
By a comparison principle for such inequalities we obtain that T{t) < T{t) and V{t) < V{t) (interpreted 
componentwise) for all t in the intersection of the domains where the solutions are defined. Here, {T{t), V{t)) 
is the solution to the linear system whose vector field appears in the right hand side of the above inequalities, 
hence these solutions are defined for all t > 0. But then T{t) and V{t) can be extended continuously to the 
closed interval [0, Tmax], contradicting maximality of I+. 

Inequality pB]) implies that for an arbitrary small e > 0, there holds that T{t) < To + e for all sufficiently 
large t. Now consider the behavior of the quantity T + I'T* along solutions of both system pO)l — p2|) and 
(Uni) - CSl): 

— (T + I'T*) = /(T) - I'BT* < f{T) - hl'T*, 
at 

where h := minj(/3i). By continuity of / on the compact interval [0,To + e], there exists (sufficiently large) 
a > such that 

/(T) + 6T < a, for aU T G [0, To + e] . 
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Therefore, for all sufficiently large t, there holds that 



— (T + I'T*) <a-bT - bl'T* < a - b(T + I'T*), 



and hence 



limsupT(i) + l'r*(t) < -. 

t~too 



Finally, from and ([15]) follows that 



and 



linisupy(i) < -r-^NB, 



linisupy(t) < -r-^P{n)NB, 
t^oo b 



respectively, where the limsup of a vector function is to be understood componentwise. Dissipativity now 
follows by observing that all the above bounds are independent of the initial condition. □ 



Lemma 2. For = 0, let all single strain equilibria Ei, E2 
and assume that 



E„ exist for either ((221) - (d) or (fl^ - (f75|) . 

To, (21) 
and 

f'{f') < 0, for all j = \,...,n + 1. (22) 

Then the Jacobian matrices of 1^10^ — (I or H13\i — () evaluated at any of the Ei's, i — I, . . . ,n+l (where 
En+i '■= Eq) have the following properties: J{Ei) has i—\ eigenvalues (counting multiplicities) in the open 
right half plane and 2{n + I) — i eigenvalues in the open left half plane. In particular, J{Ei) is Hurwitz. 

Proof. Note that when /i = 0, the Jacobian matrix associated to both model (|10|) — (fT2|) and (fT3|) — p5|) is 
the same and given by: 

//'(T) - k'V -k'T\ 
J= KV -B KT 

V NB 

To evaluate the Jacobian at any of the EiS it is more convenient to reorder the components of the state 
vector by means of the following permutations: 



1. For i = 1, . . . , n we use (T, T*, V) (T, T* , V^, ^1,^1, . 

2. For i = n + 1 we use (T, T* , V) (T, T^, Vi,Ti, V2,... 
Then the Jacobian matrices have the following structure: 

1. For i — 1, . . . 



J{Ei) 



■ ■ , T*_i, Vi-i,T*_^_^,Vi+i, 
,T*,Vn). 



,T*,Vn). 





* 




* 


■ *\ 





Bl . 








. 










. 





. 


■ Bl^i 





. 





. 







. 










. 




. 











where 




and B 



and therefore the eigenvalues of J{Ei) coincide with those of A\ and Bl, I ^ i. Since f'{T^) < it 
follows from lemma 3.4 in [T] that the eigenvalues of A\ are in the open left half plane. The matrices 
Bl are quasi-positive, irreducible matrices, hence by the Perron-Frobenius Theorem they have a simple 
real eigenvalue with corresponding (componentwise) positive eigenvector. 
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Notice that 



and thus by (pij) that 



tr{Bl) < 0, and det{Bl) = Pai 1 



A} 



< 0, for all I > i, 
> 0, for all I < i. 



There are i — 1 unstable _B-blocks on the diagonal of J{Ei 
eigenvalue to J{Ei). 



each of which contributes one positive 



2. For i 



( 



V 







where 



andB;"+i 



-Pi 
NiPi 



'11 



Notice that by a similar argument as in the previous case, all n i3-blocks on the diagonal of J(_E„+i) 
are unstable with one positive and one negative eigenvalue. 

□ 

When /i ^ 0, the question arises as to what happens to the equilibria iJi, . . . , En+i- The previous Lemma 
allows us to apply the Implicit Function Theorem which for small positive /i establishes the existence of 
(unique) equilibria Ej{^) near each Ej. Indeed, denoting the vector field of either pO| — (|12p or — (fT5| 
by F{X,ii), we have that for all j = 1, . . . , n + 1, there holds that F{Ej,0) = 0, and under the conditions of 
the previous Lemma we also have that dF/dX{Ej,0) is invertible. It is clear that £'„_|_i(/i) = £'„_|_i(0) for 
all fJ. > 0, i.e. the disease-free equilibrium is not affected by mutations. 

The main issue is of course whether or not the remaining equilibria Ej{p), j = 1, n are non- negative. 
We study this problem next and derive results in terms of the properties of the mutation matrix Q. 

For the steady-state analysis, we will need the following Lemma which is a relevant modification of 
Theorem A. 12 (ii) in [T3|. 

Lemma 3. Let M be an irreducible square matrix with non-negative off-diagonal entries and let s{M) be 
the stability modulus of M . Suppose that there exist x,r > such that Mx + r = Q. Then the following hold: 

1. If s{M) > 0, then x^r^O; 

2. If s{M) — 0, then r — and x is a multiple of the positive eigenvector of M . 

Proof. Due to Perron-Frobenius Theorem, s{M) is the principal eigenvalue of M . It is also the principal 
eigenvalue of M' . Since M' is also irreducible and non-negative off-diagonal, there exists w > such that 
M'v = s{M)v. Equivalently, v'M = s{M)v' . Hence 

= v'{Mx + r) = s{M)v'x + v'r. 

If s{M) > 0, then both non-negative products v'x and v'r must be zero which implies x ^ r ^ 0. If 
s{M) — 0, then v'r = which implies r = 0. Hence Mx = = s{M)x so that a; is a multiple of the positive 
eigenvector of M. □ 

For convenience, we introduce the following notation. We define A{fi) := T^^ NP{fi)K and assume (by 
renumbering the strains if necessary) that the strains are numbered in such a way that the matrix v4(/i) has 
the lower block-triangular form 



(Ai{^) 



(23) 
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where each diagonal block 

A,{fi) = diag^ . . . , J-J + 

is such that Bi, i — l,...,k are irreducible with non- negative ofF-diagonal entries. The off-diagonal blocks 
Bi.j, i> j are non-negative. We note that tlie diagonal entries of A{Q) are a permutation of 

1 1 1 

We say that the strain group j is reachable from strain group i if there exists a sequence of indices i = h < 
I2 < ... < Im = j such that all matrices Bi^_^^-^j^ are nonzero. Our first result is as follows: 

Proposition 1. Let the assumptions of Lemma\^hold, then the following hold: 

1. For all suffieiently small fJ. > 0, matrix A{fi) admits n distinet positive eigenvalues given by 

1 1 1 
< < • ■ • < 



Tniti) Tn-l{n) Ti(/i) 

such that Tj(0) — Ti for i — 1, ...,n; 

2. Matrix A{fi), /i > admits a positive eigenvector [vi,V2, ■■■,Vk) if and only if - ^ is a principal 
eigenvalue of Ai{fj,), and all strain groups j >2 are reachable from strain group 1; 

3. Matrix A{^), fi > admits a non-negative eigenvector {vi, V2, Vk) for each eigenvalue such that 
is a principal eigenvalue of some diagonal block Ai{p), and s{Aj{^)) < ^ ^ for all j = i + 1, k 



such that strain group j is reachable from strain group i. The component Vj is positive (zero) if group 
j is reachable (not reachable) from strain group i. 

4- All other eigenvectors of A{^), ^ > are not sign definite. 

Proof. The first assertion follows readily because A{0) has n real distinct eigenvalues and ^(/i) is continuous 
(actually, linear) in ^. The continuity of eigenvalues with respect to jj. implies that Ti(0) = for i = 1, .... n. 

To prove the second assertion, we begin with sufficiency of the condition. Let /i > be small and suppose 
that ^ I ^ is a principal eigenvalue of and all strain groups j > 2 are reachable from strain group 

1. Since Ai{fj,) is irreducible with non- negative off-diagonal entries. Perron- Frobenius Theorem implies that 
the eigenvector vi associated with ^ ^ is positive. Since is also the principal eigenvalue of it 

follows that 

s{A,{ji)-^I) <0, j>2, 

hence {Aj{ii) — - | -^)~"^ < (see e.g. Theorem A. 12 (i) in [13 ). The remaining components V2,...,Vk of 
the eigenvector satisfy the triangular system 

= t^B2^lVl + iA2{li) - ^^I)V2, 

= /iSajWi + ^iB3^2V2 + (AsifJ.) - ^r|^/)w3, 

= tJ-Bk,iVi -\ h ^iBk,k-iVk-i + (Akiti) - 7p^I)vk, 

Solving this system recursively, we obtain 

Since the strain group 2 is reachable from strain group 1, the vector fiB2.iVi > is nonzero. Positivity of 
the matrix — (A2(/i) — ^ "j: I)^^ then implies that V2 > 0. By induction on j, it follows that Vj > for all 

j = 2, k, and hence v — (wi, V2, Vk) is a positive eigenvector. 
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To prove the converse (the necessary condition), let v = {vi,V2, ■■■,Vk) be a positive eigenvector of A{iJ,) 
and let > ] be the corresponding eigenvalue. Since (^i(m) ^ ^ ] ■. I)vi — and Vi > 0, - ] must be the 

principal eigenvalue of Ai{^) (Perron-Frobenius Thm). It remains to prove that s{Aj[^)) < - ^ for all 
j > 2. Consider j = 2, and for the sake of contradiction suppose that s{Aj{ii) — I) > 0. Since the 
eigenvalues are real and distinct for small ^ > 0, this actually imphes s{Aj{iJ,) — jr^I) > 0. Then we have 
that ^ 

holds with non-negative vectors V2 and fiB2,iVi which are both nonzero. By Lemma [3l we have f2 = 0, a 
contradiction. Hence 5(^2 (/i) — ^ ^ , /) < 0. Proceeding by induction on j, we find that s{Aj{p)-^^I) < 

for all 7 > 2. Therefore, i^rrr must be the principal eigenvalue of A(n), that is, t^ttt = 7?rr^- This proves 
the second assertion. 

To prove the third assertion, we again start with sufficient condition. Suppose that ^ ^ is a principal 
eigenvalue of some diagonal block Ai{ii), and s(ylj(/i)) < - ^ for all j = « + l,...,fc such that strain 
group J is reachable from strain group i. It follows immediately that all matrices Ai{^) — I, I < i are 
nonsingular, and thus vi = 0, I < i. The component Vi is the eigenvector of Ai{ii) — rf^] s i and it is positive 
by Perron-Frobenius Theorem. Let j = i + 1, then one of the following holds. If i + 1 is not reachable from 
i, that is, = so that 

{A,+i{fi)-^J—I)v,+i^O 

implies Ui+i = because Ai+i{ii) — ^] I is nonsingular. If z + 1 is reachable from i and s{Ai^i{ii)— I) < 
0, then 

v^+i = - ^T—-iy^fiB,+i^,v, > 0. 

By induction on j, it follows that Vj = for all j > i that are not reachable from i and Vj > for all j > i 
that are reachable from i. Hence v = (0, 0, Vi, w^+i, Vk) is a non-negative eigenvector. 

Now we prove the necessary condition of the third assertion. Let v = {vi,V2, ■■■,Vk) be a non- negative 
eigenvector of A{ii) associated with eigenvalue Tfrrr- Let Vi > be the first nonzero component of v, that 

is, V = (0, 0, fi, Ufe). Then Vi satisfies — ^ ^ I)vi — hence ^ ] must be an eigenvalue of 

AiifJ.) — ^ ^ /. Moreover, by Perron-Frobenius Theorem, it must be the principal eigenvalue and Vi > 0. 
Now consider j = i + 1 and the equation 

(Ai+i(^) - . I)vi+i + uBi+i^iV, = 0. 

The vectors i^i+i and fiBi+i,iVi are non-negative. If s(Ai+i(/^) — I) > then by Lemma[31 fxBi+i^iVi = 0. 

Since n > and 1;^ > 0, this implies -Bi+i^i = 0. Equivalently, j = i + 1 is not reachable from i. An induction 
argument concludes the proof of the third assertion. 

The final assertion of this Theorem is a simple one. Let ^ | ^ be an eigenvalue of Ai{iJ,) but not the 

principal eigenvalue and let v = (wi, W2, Wn) be the corresponding eigenvector. Since all eigenvalues of 
A{fj,) are distinct, the matrices Ai{fi) — I, I < i are nonsingular so that = 0, I < i. Then Vi must be 

an eigenvector of Ai{ii) and it cannot be sign definite due to Perron-Frobenius theorem. It follows that v is 
not sign definite. □ 

Our second result concerns the existence and the number of non-negative equilibria for the systems 
(Uni) - lEl) and US]) - (HSl) with smaU ^ > 0. 

Proposition 2. Let the assumptions of Lemma\^hold and suppose that the strains are renumbered so that 
A{fx) has the form Ii23\) . Let Ej{ii) = {Tj{^),T* {^), Vj(/U,)) denote the nontrivial equilibria of both ^10\ — <\12\ 
and ((23) - en for small fi> 0. Then 

1. Ej{fi) is positive if and only if jtj^ an eigenvalue of A{fi) with a positive eigenvector Vj. 
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2. Ej{^) is non-negative if and only if | ^ is an eigenvalue of A{^) with a non-negative eigenvector Vj . 

3. Ejifi) i R^"+i if and only if is an eigenvalue of A{fi) with eigenvector Vj which is not sign- 
definite. 

Proof. We will prove the Proposition only for system pll)l — (the proof for (fT5|) — (fTS]) is similar). Observe 
that the equilibrium relation following from p2)) . can be expressed as T*{ijl) = (NB)^^TVj{fi). Hence, the 
signs of the corresponding components of T*{^) and Vj(p) are the same. Substituting the above expression 
into pTjl and (fT^ . we find that Vj{n) must satisfy 

[^-^NP{^l)K - -J—i]v^{^,) = [A{fi) - -J—i]v^{^,) = 0. 

Tjip) Tj[^x) 

Thus for each nontrivial equilibrium Ej{ij), the quantity ^ ] must be an eigenvalue of A{^) and must 

be a multiple of the corresponding eigenvector Vj. If Vj is not sign definite, it follows that Ej{^) ^ K^"^^. 
For all Vj > 0, the components of Ej{^) are uniquely determined via 

-lj[fJ,)K Vj 

Hence Ej (fi) is positive (non-negative) if and only if Vj is positive (non- negative) . □ 

An immediate corollary to Propositions [T] and [2] is that if the mutation matrix Q is irreducible, then A{fi) 
is irreducible and systems (|10p — p2p and (|13p — (jlSp with small > admit a unique positive equilibrium 
Ei{fi) and no other nontrivial non-negative equilibria. If the mutation matrix Q is reducible, then positive 
equilibrium exists if and only the fittest strain (with lowest value Ti ) belongs to strain group 1 and all other 
strain groups are reachable from group 1, meaning that the fittest strain can eventually mutate into any 
other strain. In addition, nontrivial non-negative equilibria which are not positive are feasible for > only 
if Q is reducible. Specifically, if the strains can be numbered according to (j23p . then at most k nontrivial 
non-negative equilibria exist. One extreme case is when the fittest strain belongs to group k, in which no 
positive and only one non- negative equilibrium exist. The opposite extreme case is A; = n where A{^) is 
lower-triangular, the diagonal entries of A(fi) are arranged in decreasing order, and for any pair i < j, 
strain j is reachable from strain i. In this case, there is a single positive equilibrium and n — 1 non-negative 
equilibria. 

4.1 On uniform strong repellers 

Inspired by Thieme [15j . we make the following definition. 
Definition 1. Consider a system 

X = F{x) (24) 

on a compact forward invariant set K C M™ with a continuous flow (j){t,x). Let Kq C K be a closed 
forward invariant subset of K . Let d{x,A) denote the distance from a point x to the set A. We say that 
Kq is a uniform strong repeller in K if there exists a 5 > Q such that for all solutions (j){t,x) G K\Ko, 
\immft^ood{(t){t,x),Ko) > S. 

Theorem 3. Let H : K be a continuously differentiable function such that Ii{x) ~ Q if and only if 

X G K{). Suppose there exists a lower semi-continuous function ip : K ^ R such that 

5 = ^, Vx e K\Ko. (25) 

Suppose that the following condition holds 

(H) yxeKo, 3T>0: {^p{(l){T,x)))>0. 
Then Kq is a uniform strong repellor in K . 
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Proof. Step 1. Note that by lower semi-continuity of ip and continuity of (j), for every p G Kq we can find an 
open set Up containing and a lower semi-continuous map Tp : Up ^ (0, -|-oo) so that for every q G Up, (H) 
holds with X — q and T = Tp{q). Choose for every p G Kq a non-empty open set Vp with Vp <zUp. Then by 
lower semi-continuity of each map Tp and compactness of Vp, 

inf Tp{q) > 

qeVp 

is achieved in Vp. Since Up^Ko^p is an open cover of Kq, we may choose a finite open subcover Ui^i^,,,^„Vp-. 
Let Ti = infggy^ Tp.{q) > and set 

T = min Ti > 0. 

Note that for every p G /^o, there is a T > r so that (H) holds with x = p. That is, r is a uniform (in A'o) 
lower bound for T's for which (H) holds. 
Step 2. Let h > he given. Define 

Uh = {x e K \ 3T > T : {ijj{(j){T, x))) > h} (26) 

We claim that Uh is open. 

Fix z G Ufi. Then there is some T > r so that 

e:= (7A(0(r,z)))-/i>O. 

Then by continuity of (j) and lower semi-continuity of ip (and therefore uniform lower semi-continuity of on 
compact sets), it follows that there exists an open set Wz containing z such that for all z' G Wz holds that 

i;i(t>it,z')>ij{cPit,z))^e, ViG[0,T]. (27) 

Now since 

{^iq^iT,z))} =e + h, 

it follows from that for all z' e Wz'. 

(V'(</>(T,z'))) 

and thus that Wz C Uh, establishing our claim. 
Step 3. Define Th ■ Uh ^ [t,+oo) as 

Th{z) = inf{T > T I {'^{4>{T,z))) > h}. 

We claim that Th is upper semi-continuous. 

Fix z G Uh and let e' > be given. Then there is some T > t so that 

(7A(0(r,z))) >/i, 

so that 

T<Th{z) + e' (28) 
By the argument in Step 2, there is some open set Wz containing z, such that for all z' G Wz holds that: 

{i;{<P{T,z')))>h, 

and thus that for all z' G Wz'. 

Th{z') < T (29) 

Our claim follows by combining (|28|) and (|29p. 

Step 4- The nested family {Uh}h>Q is decreasing (under set inclusion), and forms an open cover of Kq. 
Hence, there is some h so that Uj^ covers Kq. Since K :— K \ Uh is compact, and 11 is continuous, 11 attains 
its minimal value m > on K. Choose p G (0, m) and define: 

ip = {zeK\n{z)e{0,p]}. 

Then IpCUh. 
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Step 5. We claim that every forward solution starting in Jp, eventually leaves Ip, that is: 

Vze/p,3i^>0: (t){t^,z) ^ Ip. 
By contradiction, if (^(t, z) £ Ip for all t > 0, then z) G [7;^ for all t > 0, and thus: 



1 



t+Tt 



> T : — / V'('/'(s,2))ds > /i. 
Then integrating equation ([25| from t to i + T yields that: 

and thus that 

n(0(i + Tt, z)) > e''^m{(f>{t, z)). (30) 

Set to = and tk = tk-i + Tt^_-^ for k — 1,2,.... Since each Tt^ > r > it follows that tk oo. Then by 
((30)) and since t^ > t for all fc, we have that: 

n((/.(tfe,2)) > e''^*'=-in(0(ife_i,z)) > e^m(z), 

so that Il{(j){tk, z)) — > c» as A; — > cxo. This contradicts boundedness of 11 on the compact set K. 
Step 6. Let 

ip^Ip\JKo. 

We will show that there is some q G (0,p) so that forward solutions starting outside Ip, never reach Iq, that 
is: 

3q G (0,p) : z^ip^ z) i Ig, Mt > 0. 

Consider a forward solution (j){t, z) with z ^ Ip. If (j){t, z) ^ /p for all t > 0, then we are done since Iq C Ip, 
so let us assume that for some t^ > 0, holds that (/>(tz, z) ^ Ip. Denote the first time this happens by ^o- 

to = iam{t > I (j){t, z) G ip}. 

Set z* = (j>{to, z) and note that n(z*) — p. Denote infzg^^ i^iz) by m'. If m' > 0, then implies that 
n((/)(t, z*)) > n(z*) — p for all i > 0, so that we're done. If on the the other hand m' < 0, we first define 

T' = ma_xT,,(z)(>r>0). 

Notice that this maximum is indeed achieved on the compact set Ip, since Th is upper semi-continuous. Now 
we define 

and notice that q is independent of the chosen solution z{t). We will show that for this choice of q, our claim 
is established. 
We have that: 

1 /■* 

VtG(0,T): -/ ^P{(j){s,z*))ds>m', 
i Jq 

and thus by that 

ViG(0,T): n(<?!>(t,z*)) > n(z*)e™'* > g, (31) 

which implies that during the time interval {Q,T), the solution 4){t,z*) has not reached Iq. On the other 
hand, during that same time interval (0, T), the solution (j){t, z*) must have left Ip. If this were not the case, 
then by the argument in Step 5, there would be some T* G [r, T) so that 

n{^{T*,z*))>T\{z*)e^'^ >p, 

and thus that (f>(T*,z*) ^ Ip, a contradiction to our assumption. This process can be repeated iteratively 
and leads to the conclusion that the forward solution (j){t, z) which did not start in Ip, will never reach Iq. 

So far we have shown that for any solution (j){t,x) ^ Kq, inequality Il{ip{t,x)) > q> for all sufficiently 
large t. The sets Ko — n^^({0}) and n^^([q, +00)) fl K are compact and disjoint. Therefore, there exists 
6 > such that d{(j){t, x), Kq) > 6 for all x ^ Kq and all sufficiently large t. □ 
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4.2 Global stability for small yU > 

The following Lemmas will be used to prove global stability of the positive equilibrium for small /i > 0. 
Lemma 4. Let a : R" R" be continuous and let b G int(R!^). Let / : M™ x (M!^\{0}) R be defined as 

a'{x)y 

Then 

liminf /(x,?;)= min ^^j^, (32) 

furthermore, if we define f{x,0) = minjgj]^ „} , then f{x,y) becomes a lower semi- continuous function 
on R™ X R" whose restriction on R"* x {0} is continuous. 

Proof. Extending the function f{x,y) by defining /(a:o,0) = \\m\a.ix^x^^y^Q+ f{x,y) clearly produces a 
lower semi-continuous function. Furthermore, since a{x) is continuous, the function miujg^j^^ „} is 
continuous as well. So it remains to show that (15^ holds. 

Without loss of generality, we may assume that min,jg{]^ „} o-^i^o) _ ^ Setting x = xq and y2 = 

2/3 = •■• = 2/n = and letting 2/1 ^ 0+, we find that /(xq, 7/1,0, ...,0) -> ^il^ni. Hence, liminf ^^2,3,^^0+ f{x,y) < 



ai{xo) 
bi ■ 



We also observe that as long as y G R"\{0}, the value 

a'{x)y _ ^ ai{x) b^yi 



^'y ibi hyi^ 1- bnyn 

2^^, i = 1, ...,n. By continuity of cl(x), for any £ > there 
exists 6 > such that Vz S {1, n} and Va; e i?5(a;o), we have ai{x) > ai{xo)—ebi. Hence, for all x G Bs{xo) 
and for all y G R+\{0}, f{x,y) > ^^f^ - e. We have established that 

— ^ — ^> limmf f{x,yj> — e. 

bi x^xo,y^O+ bi 

Since e > is arbitrary, ([5^ follows. □ 
Lemma 5. Suppose that (0) holds. Then there exist ri,fj,o> such that 

liminf >rj>0 

t — >oo 

for any e [0, ^0] and for any solution of (f7Q|) - (HH flTid ([III - (f75|) wit/i > 0. 

Proof. We will prove the claim for system (fTO)l — (fT2|) (the proof for (fT3l) — (|15p is similar). The proof consists 
of two parts. We first show that there exist 770, /^o > such that liminft^oo l'{T*{t) + V{t)) > rjo > for all 
solutions with T*{t), V{t) ^ 0. 

We choose n positive numbers Ni so that < Ni < Ni for all 1 < i < n. This is possible since we 
assume Ti — -j^^ < To. Let v = {N, 1). It follows that 

,/-B KTo 

" [nb T 



= 61 (iVi - TVi ) , . . . , 5„ (7V„ - iV„) , /ci ToiVi - 71 , . . . , fc„To7V„ - 7, 



is a positive vector. By continuity, there exists a /io > such that 

«'M(r,M), where M{T,^^)■.^^^l Pil^)KT\^ 

is a positive vector for all ^ E [0, ^o]- 
Consider a system 



T = /(T) - fc'\/r, T e M+ (33) 

T* = P{fi)KVT - BT*, T* eWl, (34) 

F = iVBT* - Ft/, F e R+, (35) 

/i = 0, /i e [0,/io]- (36) 
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Let K' be the forward invariant compact set for PTI| — (fT^ established in Lemma[T]and define K = K'x[0, /iq] . 
It is clear that K is compact and forward invariant under (|33p — ([55]) The set 

:= ([0, To] X {0} X {0} x [0, ^o]) n K 

is clearly a compact forward invariant subset of K. 

Let Il{T*,V) :— v'{T*,V). The function H is clearly smooth, zero on Kq, and positive on K\Kq. 
Furthermore, 

n v'M{T,fi){T*,V) 

n i;'(r*,y) 

is lower semi-continuous on K by Lemma |4] once we define the value of on Kq as 

ip(l,li)— mm . 

i— l,...,n 'f;^ 

We note that the function 4'{T,n) is continuous in (T, /i). Since all solutions of ([55)1 — ([5S|) in Jfo have 
the property that limt^ooT{t) — Tq, it implies that 4'{T{t),fJ,) > for all sufficiently large t. Hence by 
Theorem 131 the set Kq is a uniform strong repeller in K. If we use the L^-norni of (T*, V) as the distance 
function to Kq, we find that there exists an 770 > such that 

liminf l'(r* + y) > r]Q 

for all solutions of ([33]) - §^ in K\Kq. 

To complete the proof, we need to show that there exists 77 > such that liminft_>oo ^'V{t) > 77 > for 
aU solutions with I'Vit) > 0. Observe that l'V{t) > implies that l'T*{t) > 0. Hence by the result of 
part one, we have that hminff^oo + V{t)) > r/o > 0, or equivalently, VT*{t) > 770/2 - VV(t) for all 

sufficiently large t. We substitute this inequality into ([72]) and find that 



I'V >Aq(^^- l'V{t)^ - Ail'V{t), Aq := mm{N,l3,) > 0, A^ max(7,) > 
holds for large t. It follows immediately that 

□ 

Lemma 6. Let 

, , z 1 

a[x, y,z) := X + y 3z3 . 

xy 

Then for any zq, M > 0, there exists (5 > such that a{x, y, z) > M for all < x < S, all y > 0, and all 

Z > Zq. 

Proof. Observe that the minimum of the function cr(x, •, z) on the set y G (0, +00) is achieved ai y = z/x. 
Hence for all 7/ > 0, it holds that 

f{x, y, z) > f{x, z)^x + - 3z3 . 

Let Zq > and define 

d :— 



{M+3ziy 

Then for all < a; < (5, all 7/ > 0, and all z > zo, it holds that 



fix, V, z) > 2^ - 3z^ ^ -| - 3z-M > zl (J^-3z,--]^ M. 



□ 
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Theorem 4. Let K he the absorbing compact set established in Lemma{^ and let 

u = {{T,T*,v) e R+"+^|r,ri*,yi > o}. 

Suppose that (C) holds with Ti instead ofT. Then there exist /ii > and a compact set Kg C U such that 
for any /i e [0, fii] and for any solution of (| — l\12^ or (jJj'P — (|_?5p in U , there exists a t^ > Q such that 
{T{t),T*{t),V{t)) e Ks for all t > to. 

Proof. Both for system (fTO|) — (fT2| and (fT3l) — p^ . the proof will be based on the same Lyapunov function 



W 



T / fl 

1 - — I dr + 



that we used to show competitive exclusion with /i = 0. 
Case 1: System (fTO|) - 

Computing W for system ([T0| — (fT2|) . we obtain after some simplifications 



^ + 



T T^ViT^ ViT* 



^ /c,y,(r^ - Ti) + -±—^^,Y.q^^k,V,T + l^Y.T. I'^okjV.T. 

i=2 1 j=l i=2 j=l 



Recombining the terms, we further obtain 



W 



(/(T)-/(T^))(1- — ) -/3iT; 



T T*ViT^ ViT* 



i=i 



i=l 3 = 1 



We note that 



E E = E ( E -^^O ^^-^^^ = 

i=l j = l j = l 



since all column sums of Q are zero. Hence, 



W 



(/(T)-/(r^)) 1 



r TlViT^ ViT* 



T* 



J=2 



We rewrite W as 



W 



(/(r)-/(fi)) 1 



T 



h (1 + qiiu) - -1 + — — ^ 

i=2 1 j=2 



Note that the last term of W is non-positive, henccfl] 

fl 



W < (/(T)-/(T^)) 1 



T 



3(l + giiM)'/' 



3(l + giiM)'/' 



+3/3ir*(i - (1 + giiAi)'/') - E ^*^'(^'' - ^')- 



1=2 



^Incidentally, if gn = 0, we obtain global stability of the boundary equilibrium Ei for all fi > 0. 
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By Lemma [SI there exist 77, /ia > such that l'V'(i) > 1] for aU /i G [0,/ia] and all sufficiently large t. Let 
a = mini>2 ki{T^ — T^) > 0, then 

n n 

J2k^V^{r -f^)>aY,V^>a{v-Vl). 

1=2 1=2 

Thus, by shifting time forward if necessary, we have the inequality 



W < {f(T)-fif')){l-^]-P^T* 



+3Pif*{l - (1 + giiA*)'/') - a?7 + 
Let /lb > be such that for all /i G [0, /if,] , 

'1 



1 + qll^J. G 



-, 1 



, 3/3ir*(l-(l + gnM)'/')-a^< 



Let //I = min(/ia, //{,) and choose sufficiently large L > so that 

3/3iT*(l - (1 + qiifiy/^) ~ar] + aVi<L 
for all solutions of (|10p — (fT2|) in /v and all /i G [0, /ii]. For any /i G [0, fii], we have that 



W < {f{T)~f{T'))[l-^]-f3^T* 



T 



Ti T*ViT ViT^ 



3(l + giiM)i/3 



at] 



aVi, 



where the first two terms are non-positive and 1 + qiil^ G l] . Inspecting the first term in W, we find 
that there exists Sq > such that 



(/(T)~/(Ti) (^1-^^ <-(L + l) 



for all T < So and all fi G [0, /ii]. Now we inspect the the second term in W. Using Lemma[6]with 

ViT: T1 ^ 1 

^^wr*' y = ^> z = i + gii/i, ^0 = 2' 



we conclude that there exists 61 > such that 



for all -y^ < 5i and all /i G [0, /ii]. Finally, there exists 82 > Q such that 



<-(L + l) 



a^i < for all Vi < S2 



and all fi G [0,/ii]. Let 

Ks = {(r,T*,i/) G ii:nc/|r > 5o,Vi > 62, t* > SiVi}. 

Consider (T, T*,V) G (iiTn U)\Ks and let /x G [0, /fi], then at least one of the following holds: 

• T < 6q, in which case W < -{L + 1) + L < -1; 

• T*/Vi < Si, in which case W < -{L + 1) + L < -1; 

• Vi < 52, in which case W < 
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Hence, for all (T, T*,V) e{Kn U)\ks and all 6 [0, fii], we have 



W < -min(l,^) < 0. 

We postpone the rest of the proof until we have showed that a similar inequality holds for system ^3]) — ([15] 
Case 2: System ^ - (fTS]) . 

Computing W for system ([T3| — (fT5|) . we obtain after some simplifications 



W = {f{T)-f{f'))(l- 



T 



/3if* 



Vi 



T T*ViT^ ViT* 

N. 



i=2 \ ^ / 

Note that the /i dependent terms can be rearranged as follows: 



i=2 j = l 



N. 



N,, 



In the above the first term is non-positive, and the second term can be re-written as follows: 



for suitable a; > 0, and the third term will be absorbed in the square bracket [ ] term in W. We find that 



W < (/(T)-/(Ti)) 1 



T 



- I3if* 



ViT* 



H -., + (1 + Qllli)- 



3(l + giiM)'/' 



-f3/3if *(1 - (1 + giiA^)'/') -f mE - E - ^')- 



i=2 



By Lemma [51 there exist r/, /ia > such that l'V{t) > rj for all fi G [0, fia] and all sufficiently large t. Let 
a = minj>2 h{f' ~ f^) > 0, then 

n n 

1=2 1=2 

Thus, by shifting time forward if necessary, we have the inequality 



W < [f{T)~f{T'))\l~I^]-P^Tl 



ViT* 



3(l + giiAi)i/3 



+3Pif*{l - (1 + giiAi)'/^) + mE - + "^1- 



Since solutions are in the compact set K for sufficiently large times, there is some /z^ > such that 

n 

M5^a.I^*<^, VMG[0,/i'J, 



and therefore 



< (/(T) - /(Ti)) ( 1 - — ) - /^iTi' 



1 --, + (1 + giiAi)- 



3(l + giiM)i/3 



-f3/3if;(l - (1 + gii/.)i/3) -^+aV^. 
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Let /^f, > be such that for all /i G [0, yUf,], 



1 + qii^i e 



1 



Let /ii = min(/ia, /ijj, /ih) and choose sufficiently large L > so that 

3/3iT*(l - (1 + qiifi)^^^) ~arj + aVi<L 
for all solutions of (|13p — in iiT and all /i e [0, /ii]. For any fi e [0, fii], we have that 



< {f{T)~f{T'))[l-^]-f3iT* 



T 



1 + (1 + giiAi)- 



3(1 + 911^)1/3 



where the first two terms are non-positive and 1 + gn/i G [5, l]- Inspecting the ffi'st term in we find 
that there exists 5o > such that 



(/(T) - /(Ti) 1 



<-(L + l) 



for all T < (5o and all /i G [0, /ii]. Inspecting the second term in W , we use Lemma [S] with 



a; = (1 + gii^) 



ViT* 



and conclude that there exists 61 > such that 



-AT* 



^ . +(1 + 9iim):^ -3(l + gii/i)'/' 



<-(i + l) 



for aU ^ < (5i and aU fi G [0, ^1]. Finally, there exists 62 > such that + aVi < for aU Vi < 82 
and all /i G [0,/ii]. Let 

Ks = {{T, T*,V) eKnU\T> do,Vi > S2, T* > SiVi}. 
Consider (T, T*,V) G {K D U)\Ks and let /i G [0, ^1], then at least one of the following holds: 

• T <Sq, in which case W < -{L + 1) + L < -1; 

• T^/Vi < Si, in which case W < -{L + 1) + L < -1; 

• Vi < 82, in which case W < 

Hence, for all (T, T*,V) e{Kn U)\ks and aU ^ G [0, m], we have 

W < -min(l,^) < 0. 

8 

The remainder of the proof is the same for both of the above two cases and presented next. 
The non-negative function W{T, T*, V, /i) is continuous and bounded from above on the set Kg x [0, /ii] 
because T,T*,Vi are bounded away from zero. Hence it attains a finite positive maximum 

w:= max W(T,T* ,V, fi) > 0. 

KsxlO,t^i] 



Define a new set 



Kg = {{T,T*,V) eKnU\W{T,T*,V,fi) < u;,V/i G [0,fii]}. 



By construction, we have that Ks C Ks C K O U. The continuity of W implies that Ks is closed, and 
therefore compact in U. 
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It remains to show that all solutions of ^TU\i — (|12p in U enter and remain in Kg for all sufficiently large 
times. Since K nU is an absorbing set for all /i > (Lemma [T|), without loss of generality we need to prove 
this for all solutions in K nU. 

Let $(<) = {T{t),T*{t),V{t)) e KnU he a solution of ((TOl) - (HH) for some fixed fi e [0,^i]. Observe 
that in the set {K n U)\ks, the inequality W < - min(l, ^) < holds. Since W > 0, there exists to > 
such that $(to) G Kg C Kg. We will show that € Kg for all t > to. For the sake of contradiction, let us 
suppose that there exists ti > to such that ^{ti) ^ Kg. Then there exists t2 £ [to,ti) such that $(^2) G Kg 
and ^ for all t £ (^2,^1]- On the one hand, we have that 

w{<^>{t2),^i)<w<wmtl),^i) 

by definition of Kg. On the other hand, for all t G (^2,^1], we have $(t) ^ Kg and consequently <i>(i) ^ Kg 
so that ^T4^(<f>(t), ^) = < 0. This contradiction shows that G for all t > to and concludes the 
proof of the Theorem. □ 

Theorem 5. Let the assumptions of Lemma\^hold, let U be the set from Theorem^ and define 

U' = {(T, T\ V) e R^"+i| T* + > 0} D [/. 
Then there exist /io > and a continuous map E : [0, /io] ^ such that 

1. E(0) — El (where Ei is the same as in Lemma\^, and E(fj,) is an equilibrium of (jtOp — (jj^p or of 
((23) - en /or all n e [0,/io]; 

2. For each fi G [0, /lo] , is a globally asymptotically stable equilibrium of l\10^ — l\12^ or of (|_? j| — (|_?5p 
m [/'. 

Proof. To prove the first assertion, we begin by noting that for yu = 0, -Bi is a stable hyperbolic equilibrium 
of ITOD - dig) or of - (dH) by Lemma H Since the vector field of ITUl) - (HI]) and US]) - dH]) is linear 
in /i, by the Implicit Function Theorem there exist h > and a continuous map E : {—h, h) M^"+^ such 
that E{p) is an equihbrium of HO]) - dH]) or ^ - ^ for all fj. G (-/i, /i). The fact that E{p.) G C/ for all 
/i G [0, /i) follows from Proposition [2] and the fact that fi < Ti, i > 2. Note that for fi > 0, E{ii) may be 
positive (if Q is irreducible) or non- negative (if Q is reducible). Nevertheless, in both cases, /i > implies 
E(pl) G U. 

The proof of the second assertion is based on the result of Smith and Waltman (Corollary 2.3 in [14 ). 
We have already established the fact that E{Q) is a stable hyperbolic equilibrium of (fTU]) — (fT^ or (fO)) — ([T5|) . 
By Theorem[2l -E(O) is globally asymptotically stable in U' for ji — 0. In addition, by Theorem |4] there exist 
/^o > and a compact set Kg C U such that for each ^ G [0,^o]; and each solution (r(t), r*(t), y(t)) of 
dTO]) - (ini) or - (HH) in [/, there exists to > such that {T{t),T*{t),V{t)) G Kg for all t > to. Hence, 
the condition (HI) of Corollary 2.3 in [T3] holds. The Proposition 2.3 itself then implies the global stability 
of E{fi) in U for all sufhciently small /i > 0. Finally, solutions of (dU]) - dH]) or dHI - starting in U' 
enter U instantaneously, hence global stability of E{ii) in U' follows as well. □ 



Appendix: Inclusion of loss of virus in the model 
Single-strain 

When taking the loss of the virus particle upon infection into account, model dH) becomes 

f = /(T) - kVT 
f* = kVT - f3T* 

V = N(3T* -jV - kVT, (37) 

We still assume that the growth rate of the healthy cell population is given by dH), hence Eq — {To, 0, 0) is 
still an equilibrium of (j37p . A second, positive equilibrium may exist if the following quantities are positive: 

rfi T rp^ f (--^) / (-^) (38) 



k{N-l)' (3 ' kT 
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Note that this is the case iff iV > 1 and / ^ fc(jv-i) ^ > 0; or equivalently by ([2]) that T = ^.^^^^^^-^ < To- In 
terms of the basic reproduction number 

0. k{N-l) To 

existence of a positive equihbrium is therefore equivalent to T?." > 1. Assuming that TZ'^ > 1, we wiU stiU 
denote this disease steady state by = (f , T*,V). We introduce the foUowing condition. 



(C) /'(c) + -f{f) < 0, for aU c e [0,To]. 

7 

Note that this condition is satisfied when /(T) is a decreasing function with sufficiently large negative 
derivative. 



Theorem 6. Let (C) hold. Then the equilibrium E is globally asymptotically stable for l\37l with respect to 
initial conditions satisfying T*(0) + V(0) > 0. 

Proof. Consider the following function on int(Mi^): 



T 



W=iN-l) [l--]dT + N 



T* 



1 - f 1 dr. 



Then 



W = {N-l){f{T)-kVT)(^-^^+N{kVT-pT*)(^-^^+{NI3T*--iV-kVT)(^-^ 
= (iV - 1)/(T) (l - I) - NkVT^^ + NPT* - NI3T*^ + 7? + kVT 



= {N- l)(/(r) - fit)) {l-T^+{N- 1)/(T) (1 - ^) + NPT* 

= {N- l)(/(r) - /(T)) (^1-I.^+(N~ 1)0* (1 - ^) + Nf^T* 

= (TV - 1)(/(T) - /(T)) (1 - ^) + (^-2 + ^ + 1^ + 7V/3f * 

= [{N - 1)(/(T) - /(r))T + (3f*{T ~ T)\ ^^^^ + Nlif* 



2 - 





T*V' 


YTT* 




V _ 


VTT* 


T*V' 


VTT* 


T*V 


VTT* 


T*V 


VTT* 


f*V 


T*V 


T' 




f*V 


~ T 





- (3T* + f3T*'^ 
(3f* +pf*'^ 



where we used (1381) repeatedly; in particular in the second, third and fourth equation. By the mean value 
theorem there is some c G (T, T) or (T, T) such that 



hence using (|38p once more 

W ={N -I) 



f{c) + -f{f) 
7 



/(T)-/(r) = /'(c)(r-r), 
(T - ff 



T 



■N/3T* 



VTT* T*V T 



VTT* T*V T 



The first term is non-positive by (C) and because we can assume that T < To by dissipativity (see Lemma 
[7] later). The second term is non-positive as well since the geometric mean of 3 non- negative numbers is not 
larger than the arithmetic mean of those numbers. We conclude that < in int(Mi^), hence local stability 
of E follows. Notice that W equals zero if and only if both the first term and the second term are zero, This 
happens at points where: 

f f*V 

— = 1 and ^ = 1. 

T T*V 

Then LaSalle's Invariance Principle implies that all bounded solutions (and as before, solutions are easily 
shown to be bounded, see also Lemma [7] later) in int(R'^) converge to the largest invariant set in 

M = {{T,T*,V) e int(R3 ) | ^ = 1, |^ = 1}- 
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It is clear that the largest invariant set in M is the singleton {E}. Finally, note that forward solutions 
starting on the boundary of R"^ with either ri(0) or Vi(0) positive, enter int(M5j_) instantaneously. This 
concludes the proof. □ 

Competitive exclusion 

Now we modify the multi-strain model ([4]) — ([6]) to 



f = f{T) - kVT, T e M+ 
t* = KVT-BT*, T* eRl 
V = NBT* - TV - KVT, V € 



(39) 
(40) 
(41) 



where k = (fci, fc„), K = diag(fci, ...,fc„), B = diag(/?i, ...,/?„), N = diag(iVi, ...,iV„), andT = diag(7i, ...,7„)- 

itate by itself, that 
, En- Assume that 



Suppose that each strain is capable to persist at steady state by itself, that is, Ni > 1 and Ti = < To 

and denote the corresponding equilibria also by Ei, 

0<Ti<T2<...<T„<To. (42) 

In addition, suppose that (C) holds with T = Ti. Then we have the following. 

Theorem 7. The single strain equilibrium Ei is globally asymptotically stable for l\39\i — with respect 
to initial conditions satisfying T]*(0) + Vi{0) > 0. 



Proof Consider the function W defined on U := {{T,T*,V) G 



p2n+l 



T,Tl,Vi > 0} as 



W = (iVi - 1) 



1 - — ) dr + 



TT 



Vi 



" AT - 1 



Computing W, we find that 
W = {Ni-l 



f'ic) + -/(Ti) 

7 



(T - fif 



J2 ( -hViT-Ti, 

i=2 



T 

(Ni - 1) 



NiPiT* 



ViTT* T*Vi Ti 



ViTiT* T*Vi T 
{N.hVT - NAT* + N,(3,T: - 7, - hVT)^ . 



After simplifications, we have 

W = (A^i - 1) 



/'(c) + -fin) 

7 



(T - fi)2 



T 



ViTT^ T^Vi Ti 
WT\f* ~ T*Vi ~ Y 



-{Ni~l)Y,hV,{f,-fi). 



The first term is non-positive since (C) with T = Ti holds and because T < Tq by disspiativity (see Lemma 
[7] later) . The second term is non-positive is well, and so is the third by ((42|) . Thus W < which already 
implies that Ei is stable. An application of LaSalle's Invariance Principle shows that all bounded solutions 
in U (boundedness follows from Lemma [7] which is proved later) converge to the largest invariant set in 

(T,r*,...,T„*,yi,...,K)ec/| 1^ = 1, = v^^o, i>2 

which is easily shown to be the singleton {Ei}. Finally, solutions on the boundary of U with rj*(0)+Vi(0) > 
enter U instantaneously, which concludes the proof. □ 
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Adding mutations 

We modify the model (|39p — (HT]) to account for mutations. Again, we consider two alternative models 

f = f{T) - kVT, T e M+ 
f* = P{fi)KVT - BT*, T* eWl 

V = NBT* - TV - KVT, V eR^, (43) 

and 

T /(T) - kVT, T G R+ 
T* = KVT-BT*, r*GR" 

V = P{ii)NBT* - TV - KVT, V eW^, (44) 

where A:, if, B, N , T are the same as before, and P(/x) = I+nQ and Q is a stochastic matrix with non-negative 
off-diagonal entries. 

Lemma 7. Both systems \4^S^ and are dissipative, i.e. there is some compact set K such that every 
solution eventually enters K and remains in K forever after. 

Proof. The proof is similar to the proof of Lemma [T] and will be omitted. □ 

Lemma 8. For ii — 0, let all single strain equilibria Ei, E2 ■ . . , En exist for either ^J3^ or ^44h '^^^ assume 
that 

< f 2 < • • • < f " < f To, (45) 

and 

f'(f^) < 0, for all j = l,...,n+l. (46) 

Then the Jacobian matrices of or ^44h evaluated at any of the Ei's, i — 1, . . . , n+1 (where En+i '■= Eq ) 
have the following properties: J{Ei) has i — \ eigenvalues (counting multiplicities) in the open right half plane 
and 2{n + 1) — i eigenvalues in the open left half plane. In particular, J{Ei) is Hurwitz. 

Proof. The proof is similar to that of Lemma[2] The only difference is that the entries of the Jacobian matrices 
change. In particular, the (3, 1) and (3,3) entry of A\ now become —kiVi and —ji — kT^ respectively, but 
by (|¥5|) and Lemma 3.4 in [T, A\ is still Hurwitz. 

□ 

To study equilibria of systems (|43|) and (|44p . we introduce the matrix 

A{^Ji)=T-\NP{^i)-I)K, (47) 
which has non-negative off-diagonal entries for /i > and 

^,o) = aiaJM^,...,M«W)U /I ly 

V 71 In J \Ti TrJ 

Clearly, Proposition [T] holds with A(/i) given by (|T7)) . Hence, we have the following. 

Proposition 3. Let the assumptions of Lemma\^hold and suppose that the strains are renumbered so that 
A{fi) has the form I123\) . Let Ejlp,) = {Tj(fj,),T*(fi),Vj{fj,)) denote the nontrivial equilibria of both l \43l and 
l\44l for small fi > 0. Then 

1. Ej{ii) is positive if and only if ^ ^ is an eigenvalue of A{ii) with a positive eigenvector Vj . 

2. Ej{ii) is non-negative if and only if | ^ is an eigenvalue of A{fi) with a non-negative eigenvector Vj . 

3. Ejifi) i R^"+i if and only if ^ - is an eigenvalue of A{fj,) with eigenvector Vj which is not sign- 
definite. 
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Proof. We will prove the Proposition only for system (the proof for is similar). Observe that at 
equilibrium, T*{ii) = {NB)~^{T + KTj{fi))Vj{fi). Hence, the the signs of the corresponding components of 
T*{fi) and Vjin) are the same. Substituting the above expression into ((43|l . we find that Vj{iJ,) must satisfy 

[r-i(7VP(M) - I)K - -J—I]V^{^) = [A{p) - -J—I]V^{f,) = 0. 

Thus for each nontrivial equilibrium Ej{fi), the quantity ^ ] must be an eigenvalue of A{^) and Vj{p) must 

be a multiple of the corresponding eigenvector Vj. If V^- is not sign definite, it follows that Ej{^) ^ R^"^^. 
For all Vj > 0, the components of Ej{^) are uniquely determined via 



Hence Ej (fi) is positive (non-negative) if and only if Vj is positive (non- negative) . □ 



Lower bounds 

Lemma 9. Suppose that l\45^ holds. Then there exist rj, iiq > such that 

liminf lV(t) > t] > 

t — ^oo 

for any fi G [0, /io] and for any solution of l\43\j and with l'V{t) > 0. 

Proof. We will prove the claim for system fthe proof for (|44p is similar). The proof consists of two parts. 
We first show that there exist 770, Mo > such that liminft^oo l'{T*{t) + V{t)) > 770 > for all solutions 
with T*{t), V{t) 0. We choose n positive numbers Ni so that < < Ni for alll < i < This is 

possible since we assume Ti = k {N ~i) -^0 which is equivalent to Ni > ■ Let v = {N, 1). It follows 

that 



, f-B KTo 
\NB -r-KTn 



h{Ni - iVi), 6„(7V„ - iV„), kiToNi - (71 + fciTo), k„ToN„ - (7™ + knTo) 



is a positive vector. By continuity, there exists a fiQ > such that 

v'MiTo,n), where M(T,m) 1^;^;^ ^T^I^Kt) 

is a positive vector for all ^ G [0, /io]- 
Consider a system 

f = f{T) - k'VT, T G M+ (48) 
T* = P{fi)KVT ~ BT*, T* eWl, (49) 

f = iVBr*-rF-ii:Fr, fgm+, (50) 

fi = 0, /i G [0,Mo]- (51) 

Let K' be the forward invariant compact set for (|43)) established in Lemma[7]and define K ^ K' x [0, /io]. It 
is clear that K is compact and forward invariant imder — (|51l) The set i^o = ([0, To] x x x [0, /io]) Hii' 
is clearly a compact forward invariant subset of K. 

Let H(T*,V^) := v'{T*,V). The function H is clearly smooth, zero on Kf), and positive on K\Kq. 
Furthermore, 

H ^ ^ v'M{T,ti){T\V) 

H ^'(r*,^) 

is lower semi-continuous on K by Lemma 2] once we define the value of ^ on i^o as 

ip(l,li)— mm . 

i=l,...,rt i;,- 
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We note that the function tp{T,fi) is continuous in (T, ^). Since all solutions of (gS]) — ([?T|) in Kq have 
the property that limj^oo T(t) = Tq, it impHes that ip{T{t),ii) > for all sufficiently large t. Hence by 
Theorem [21 the set Kq is a uniform strong repellor in K. If we use the L^-norm of (T*, V) as the distance 
function to Kq, we find that there exists an 770 > such that 

liminf l'(r* + y) > rjQ 

for aU solutions of (gS]) - (gl]) in K\Kq. 

To complete the proof, we need to show that there exists 77 > such that liminft^oo l'^(i) > ?7 > for 
all solutions with l'V{t) > 0. Observe that l'V{t) > implies that l'T*{t) > 0. Hence by the result of 
part one, we have that hminft_oo + V{t)) > r/o > 0, or equivalently, l'T*{t) > 770/2 - l'V{t) for all 

sufficiently large t. From ([50]) . we have that 

n n n n 

I'V > N^^3^T: ~ ^(7, + hT)V, > N^P^T* - ^(7, + hTo)V,. 

i=l i=l i=l i=l 

Hence, 

I'V > Ao f ^ - l'V{t)] - AilV(t), Aq := min(A^,A) > 0, Ai := max(7, + fc.To) > 



^ 2 

holds for large t. It follows immediately that 



liminflV(0>^= .,7;^° ^ 
*^oo 2(^0 + Ai) 

□ 



Existence of absorbing compact set for small /i > 0. 

Theorem 8. Let K he the absorbing compact set established in Lemma^ and let 



u ^ {{T,T*,v) e R+"+^|r,ri*,i/i > o}. 



Suppose that there exists e > such that 



ki 



(C) /'(c) + — /(fi) < -e < 0, for all c e [0,Tq]. 

71 

Then there exist > and a compact set Kg C U such that for any /i G [0, /ii] and for any solution of 
system in U, there exists a to > such that (T(t), T*{t),V{t)) £ Kg for all t > to- 
An identical statement holds for system <\44h 



Proof, (a) We first prove the statement for system (j43|). Consider the function 



W 



(iVi - 1) 



fi 



dr + Ni 



T' 



dr- 



Vi 



dr 



" AT - 1 



i=2 



Computing W for the system (|43)) . we obtain 



W = (iVi - 1) 



/'(c) + -/(Ti) 

7 



(T-Ti 



T 



ViTiT* T*Vi T _ 

(iVi - 1) ^ k,V,{f, - Ti) + /i^ /^ X] 91,^,^,^ + - 1) E E ^vkjV.T, 

1=2 1 1 = 1 1=2 * 7 = 1 
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Recombining the terms, we find that 



W 



{N, - 1) 



/'(c) + -/(Ti) 

7 



T 



ViTT^ T*Vi Ti 



i=2 



■ 1 * j=i 



"^TViQii 



1 j=i 



where the last term is clearly non-positive. Let 



a = (A^i - 1) min h{T, - Ti) > 0, 

By Lemma ini there exist ri,^a > such that l'V{t) > for all /x G [0,/i-(j] and all sufficiently large t. Hence, 
by shifting time forward if necessary, we have the inequality 



T 



-a{ri -Vi) + ^.L- fiNiqii 



VifiT* ~ T*Vi ~ Y 



T*ViT 



T* 



which holds in K for all fi G [0, /ia]- We combine the second and the last terms to obtain 



Further, we rewrite the above inequality as 



3 - (1 + qiin) 



ViTiT* T*Vi T 



T 

-a{f] -Vi) + tiL + 3Nipif* 
Let /If, > be such that for all /i G [0, /if,]. 



1 - (1 + qiifiy/^ 



(1 + giiAi) e 1], -c^v + f^L + 3Nif3if* [l - (1 + giiM)'/' 



< 



arj 



Now we let fii = min[/ta, /if,], so that for all fi G [0, /ti] and all points in 



T 



3(1 + 911/^)1/3 _ (l + giiAi) 



ViTT* T*Vi Ti 
VifiT* ~ f*Vi ^ Y 



-— + aVi. 



Let Li = asup^ Vi. Inspecting the first term in W, we find that there exists Sq > such that 

for all T < do- Similarly, inspecting the second term in W and using Lemma [6l we find that there exists 
Si > such that 



VifiT* ~ f*Vi ~ Y 



< ~Li 
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for all < Si and all £ [0, ^i]. Finally, there exists ^2 > such that — ^ + aVi < — ^ for all Vi < S2 



Let 



Ks = {{T,T*,V)eKnU\T>Sn,Vi>S2,T*>SiVi}. 
Consider (T, T*,V) G {K Ci U)\Ks and let fi e [0, ^1], then at least one of the following holds: 



• r < (5o, in which case W < —Li 



+ ii<-^; 



• Tl/Vi < Si, in which case W < -Li - ^ + Li < -fi; 

• Vi < S2, in which case W < 

Hence, for all (T, T*,V) €{Kn U)\ks and all fj. G [0, m], we have W <0. From this point forward, 

the proof is identical to the proof of Theorem IH so it will be omitted. 

(b) Now we consider system Let W be the same as in part (a). Computing W for the system pi]) . we 
obtain 



W = {Ni- 1) 



/'(c) + -fifi) 

7 



(r - Ti) 



i=2 1 j = l 1=2 ' i=l 



ViTT* r*Vi Ti 
TVi - 1 



Recombining the terms, we find that 



W 



(Ni - 1) 



/'(c) + -/(Ti) 

7 



(T-Ti) 



VifiT* ~ f*Vi ^ Y 



i=2 



-(Ni - 1) E - Ti) + E E *.^./^.^, 

=1 ' J=l 



-Wii- 



14 



i=2 



where the last term is clearly non-positive. Let 

a = (iVi - 1) min k, {f, - Ti) > 0, 



i>2 



N -1 " 

^ = supEir^E*^-^^/3,^;>0. 



By Lemma ini there exist rj^^a > such that l'V{t) > rj for all ^ G [0,/ia] and all sufficiently large t. Hence, 
by shifting time forward if necessary, we have the inequality 



W < ^e(jVi-l) ^^ J'^^ +N1P1T* 



ViTT^ T^Vi Ti 
V{T\f* ~ f*Vi ~ Y 



-a{r] -Vi) + ^iL- ^iqii 



ViNiPiT^ 



which holds in K for all /i G [0, fia]- We combine the second and the last terms to obtain 



W < -£(JVi-l) ^^ +Ni(5if* 
-a{rj - Vi) + fiL. 
Further, we rewrite the above inequality as 

W < ^e{Ni - 1) + Nif3iT* 



ViTiT* 



T*Vi T 



ViTiT* 



T*Vi T 



MV - Vi) + iiL + iNipifl [1 - (1 + qiiiif'^ 
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Let /^f, > be such that for all /i G [0, /^f,], 

arj 



(1 + giiM) e [i 1], -ar; + + SN^Pif* [l - (1 + 911^)'^' 
Now we let fii = min[/ia, /if,], so that for all /t G [0, /^i] and all points in if, 



< 



ViTTi* T*Vi Ti 



3(1 + gnM)'/'-^Wr- (1 + 911/^) 



1/1 TiT* ' "'r*yi r 



From this point forward, the proof is identical to the proof of part (a), so it will be omitted. □ 
Theorem 9. Let the assumptions of Lemma\^hold, let U be the set from Theorem\^ and define 

U' = {(T, T*, V) G M+"+^| Ti* + > 0} D [/. 
Then there exist fiQ > and a continuous map E : [0, /io] U such that 

1. E(0) = El (where Ei is the same as in Lemma\^ , and E{fi) is an equilibrium of l\43^ or of l \44i for 
all fi G [0, ^q]; 

2. For each /i G [0,/io], E{^) is a globally asymptotically stable equilibrium of ^43^ or of \4'i\ ™ U' . 
Proof. The proof is similar to that of Theorem [51 □ 
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